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It is shown that a Dirac particle of mass m and arbitrarily small momentum will tunnel without reflection through a potential 
barrier V = U c (x) of finite range provided that the potential well V = — U c (x) supports a bound state of energy E — —rn. 
This is called a supercritical potential well. 
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It is now over 70 years since the Dirac equation was 
written down. Yet new results have been discussed in 
recent years, even in the relatively simple cases of one 
|jj and two spatial dimensions Q as well as in three di- 
mensions 0. In this note, we generalise a well known 
theorem of scattering off a one-dimensional potential well 
in the Schrodinger equation to the Dirac equation. This 
is not difficult, but the theorem has an unexpected twist. 
Since the Dirac equation covers anti-particle scattering 
as well as particle scattering, the generalisation gives two 
distinct results. One of these results implies a remarkable 
property of tunnelling through a potential barrier in the 
Dirac equation which is related to the result on barrier 
penetration found by Klein and now called the Klein 
Paradox. 

We begin by considering the scattering off a class of 
one-dimensional potential wells V(x) where V(x) = 
for |x| > a and V(x) = —U(x) < for |a;| < a where the 
piecewise continuous function U(x) > . The potential is 
also taken to be even so that V{— x) = V(x). We first seek 
to generalise to the Dirac equation the non-relativistic 
result that the reflection coefficient R(k) for scattering 
off the potential well V(x) — —Uq(x) which supports a 
zero energy resonance satisfies R(0) = 0, where k is the 
momentum of the particle. This theorem was known to 
Schiff @ and Bohm g but aproof was published only 
relatively recently by Senn [|10| and Sassoli de Bianchi 
{□J . The situation where R(k) = and the transmission 
coefficient T{k) = 1 is called a transmission resonance 
H. In non-relativisti c sy stems a zero energy resonance 
(or half-bound state) p2[ is the non-trivial limit where a 
bound state just emerges from the continuum, for exam- 
ple when a square well potential is just strong enough to 
support a second bound state. 

Following an earlier paper |l3| ] we take the gamma ma- 
trices "f x and 70 to be the Pauli matrices a x and a z re- 
spectively. Then the Dirac equation for scattering of a 
particle of energy E and momentum k by the potential 
V(x) can be written as the coupled equations 



df 

+ (E - V(x) +m)g = 

ox 



(la) 



where the Dirac spinor ip = 



Eqs. (1) have simple solutions as i — > ±oo where 
V = 0. In particular, the analogue of a zero energy res- 
onance in the Schrodinger equation is a zero momen- 
tum resonance in the Dirac equation where a particle 
of zero momentum has E = m or an anti-particle has 
E = — m.[]. It is easy to see that the solution of Eq. 
(1) for E = m and V — appropriately normalised is 
/ 2m 



is ( 2^ V As in Ref. |l3| we can now write down the 

solutions of Eq. (1) for a particle of momentum k as 

E + m 



while the solution for E = —m and V = 



±oo to obtain tp 



—ik 



while an anti- 



particle of momentum k will have tp 



ik 



-ikx 



m — E 

We now set up the usual formalism for particle scatter- 
ing by the potential V(x) in the Dirac equation. We take 
the particle as incident from the left, so the amplitude 
for reflection r(k) is denned through the spinor ip{x) as 



ip(x) = 



E + m 
—ik 



ikx 



+ r(k) 



E + rn \ _ lfc; 
ik 1 ' 



while as x — > oo 



ip(x) = t(k) 



E + m \ eikx 
—ik 



(2) 



(3) 



In the Dirac equation jlj] as for the Schrodinger equa- 
tion with symmetric potentials pd|], unitarity implies 
that 



|r| 2 + |i| 2 = l; 



Im(r*t) = 



(4) 



so R + T = 1 where the reflection and transmission coef- 
ficients are given by R — \r\ 2 , T = \t\ 2 . 



|f - (E - V(x) - 771)/ = 



(lb) 



The anti-particle is described by the hole wave function 
corresponding to the absence of the state with E = — m 
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Since the potentials we consider are even, parity is con- 
served. In our two-component approach, the transforma- 
tion of a wave function under x — > — x is given [fL3[ by 



ip'(x,t) = a z ijj(-x,t) 



(5) 



It follows that an even wave function ip has an even top 
component / and an odd bottom component g. Similarly 
for an odd wave function ip, f will be odd and g will be 
even. 

We first consider an even bound state in the potential 
V(x). As x — ► ±00, its unnormalised wave function will 
be of the form 



tp(x) 



ip{x) 



m + E 

— K 



-OO 



E 



(6a) 



(6b) 



where E 2 — m 2 — k 2 . We require the potential well 
V = —Uq(x) to just bind this bound state with arbitrarily 
small k. If this is the case, then the limit k — > exists, 
and ip( x ) becomes a continuum wave function since it is 
no longer square integrable. 

We can now compare Eqs. (6) with Eqs (2) and (3) in 
the limit k -> 0, n -> 0. We obtain 2m(l+r(0)) = 2mt(0) 



1 + r(0) = t(0) ^ 



(7) 



We have written i(0) 7^ since otherwise tp(x) would 
vanish in the limit k — > 0, k — > and we would not be 
considering a zero momentum resonance. ^ 

The theorem now follows easily just as it does in the 
Schrodinger case Combining Eq (7) with the unitar- 
ity condition Im(r*t) — of Eq (4), we get Im(r*) = 



so that r(0) and £(0) are real, 
obtain 



From 



2r 2 + 2r + 1 = 1 



1 we 



(8) 



so that r(0) = or r(0) = —1. Since t(0) 7^ we obtain 
the result r(0) = and so in terms of the reflection and 
transmission coefficients 



R(0) = T(0) = 1 



(9) 



If instead we had considered an odd bound state, an 
additional minus sign must be introduced into either 
Eq. (6a) or Eq. (6b). Eq. (7) must be modified to 
1 + r(0) = — 1(0) and the subsequent analysis and con- 
clusions remain valid. Hence just as in the Schrodinger 



2 In Ref j^] we adopt a more general approach to obtain the 
results of this letter thereby avoiding the use of i(0) 7^ : 



equation, the scattering of a particle in the Dirac equa- 
tion off a potential well V = —Uq(x) which "binds" a zero 
momentum resonance corresponds to a transmission res- 
onance with zero reflection. 

We now increase the strength of the potential well from 
Uq(x) to U c {x) so that V = —U c (x) supports a bound 
state of energy E = —to. This is called a supercriti- 
cal potential and is associated with spontaneous positron 
production Jig ), |l^ ]. We can redo the analysis exactly 
as before by defining amplitudes r_ , i_ for the reflection 
and transmission of an anti-particle of momentum k in- 
cident from the left on a potential well V(x): so in place 
of Eq. (2) we have as x — > —00 



ip{x) = 



ik 
m — E 



—ikx 



while as x — > 00, we have 



ip(x) = t-(k) 



As before unitarity gives 



+ r_(fc) 



m — E 



—ik 
m — E 



li-l =1; 



Im{r*_t-) = 



(10) 



(11) 



(12) 



and the near-supercritical even bound state for x — > ±00 
is now 



if)(x) 



ip(x) 



— K 

m — E 



K 

m — E 



00 



(13a) 



(13b) 



Note again that for the odd bound state we must drop 
the minus sign in Eq. (13b). 

Repeating the analysis of Eqs (7-9) we find in the limit 
k — > 0, k — > when the antiparticle is incident on the po- 
tential well V = —U c (x) with arbitrarily small momen- 
tum that 



#_(0) = T_(0) = 1 



(14) 



where R- = |r_ \ , T_ = |t_ | 

So we see that in the Dirac equation there are two ana- 
logues of the Schrodinger result: one for zero momentum 
particles incident on a potential well which supports a 
zero momentum resonance and one for zero momentum 
particles incident on a supercritical potential well. 

We now can obtain our main result. The Dirac equa- 
tion (1) is invariant under charge conjugation: that is to 
say under the transformation 



E^-E V^-V f^g g^f 



(15) 



From Eq (14) we know that an antiparticle of energy 
E = — V m 2 + k 2 incident on the supercritical potential 
well V c (x) = —U c (x) will satisfy T_(0) = 1, that is to say 
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at arbitrarily small momentum it will have a vanishingly 
small reflection coefficient. Eq. (14) then shows that if 
we replace the antiparticle of energy E = —\Jm 2 + k 2 
incident on the supercritical potential well by a particle 
of energy E = y/m 2 + k 2 incident on the corresponding 
potential barrier V(x) — +U c (x) the particle will still 
have a transmission resonance at zero momentum, even 
though now the potential well has been replaced by a 
potential barrier. We thus obtain the theorem that where 
an even potential well of finite range is strong enough to 
contain a supercritical state, then a particle of arbitrar- 
ily small momentum will be able to tunnel right through 
the potential barrier created by inverting the well with- 
out reflection. This result was noticed a few years ago 
for the particular case of square barriers (l3| , and one 
of us (PK) has shown numerically that it was also true 
for Gaussian and Saxon- Woods potential barriers jl8). 
In the Appendix we show the behaviour of the upper 
and lower components of the wave function for particle 
scattering at zero momentum by a square and Gaussian 
barrier when the corresponding potential wells are super- 
critical. 

Conclusions 

In his original work Klein |7j discovered that a Dirac 
particle could tunnel through an arbitrarily high poten- 
tial. The generic phenomenon whereby fermions can tun- 
nel through barriers without exponential suppression we 
have called "Klein Tunnelling" The result of this 

letter shows that Klein tunnelling is a general feature of 
the Dirac equation: any potential well strong enough to 
support a supercritical state when inverted becomes a 
potential barrier which a fermion of arbitrarily low mo- 
mentum can tunnel through without reflection. We do 
not claim here that any transmission resonance at zero 
momentum must correspond to supercriticality, only that 
supercriticality leads to a transmission resonance through 
a potential barrier at zero momentum. In another paper 
fil| we shall consider the question of the conditions on 
a potential for it to possess a zero momentum transmis- 
sion resonance more generally. In three dimensions Hall 
and one of us (ND) [M have recently demonstrated that 
maximal Klein tunnelling is also associated with super- 
criticality. 

The potential step that Klein considered has patho- 
logical properties [Q. Nevertheless our result confirms 
that according to the Dirac equation a particle of low mo- 
mentum can tunnel through an arbitrarily high smooth 
potential of finite range. The reason is straightforward: 
hole states can propagate under the potential barrier. 
In terms of the particle kinetic energy T under the bar- 
rier T=E — V — m = —to — \J m 2 + q 2 where q is the 
momentum of the hole so if T < —2m, hole states can 
propagate without exponential suppression. T < —2m 
thus corresponds to penetrating under the barrier to dis- 
tances |a;| < \xk\ where V(xk) = E + m > 2m |20|. 

Appendix 

We illustrate the result above for the special cases of 



(i) a square barrier and (ii) a Gaussian barrier. First 
consider the square well potential V = —U(x) where 
U(x) = U for \x\ < a and U(x) = for \x\ > a. Then an 
unnormalised even wave function inside the well [fl3f has 
the form 



ip{x) = 



(E + U + m) cospx 



p sin px 



\x\<a (16) 



where the internal momentum p is given by (E + U) 2 = 
m 2 +p 2 . For supercriticality where E = — m we require 
the phase condition pa — Nit/ 2 where N is an integer 
tn|. The first supercritical state is thus given by p — 
p c = n/2a and correspondingly the critical potential is 
V(x) — — U c \x\ < a, where U c — m + yj m 2 + tt 2 /4a 2 
At supercriticality the wave function is thus given by 



ip(x) = —2m 





2m 



x < 



bcos(irx/2a) 
sin(TTx/2a) 



\x\ < a 



ip(x) = - 





2m 



x > a 



(17a) 



(17b) 



(17c) 



where b = 2aU c /Tt. 

Now consider a particle of arbitrarily small momen- 
tum incident on the square barrier V{x) = U c \x\ < 
a; V(x) = \x\> a. Eq (15) shows that the wave func- 
tion is obtained by interchanging the top and bottom 
components of Eq (17) thereby giving the transmission 
resonance 



ip(x) 



-2m 



2m 




x < 



sin(7ra;/2a) 
6cos(7rx/2a) 



(18a) 



\x\ < a (18b) 



2m 




x > a 



(18c) 



and the components of the wave function are shown in 
Fig.l. 

One of us (PK) has also solved the Dirac equation nu- 
merically for a Gaussian potential well and barrier where 
U(x) = U exp(—x 2 /a 2 ) Q. In Fig. 2 we show the com- 
ponents of the wave function for a particle of arbitrar- 
ily small momentum incident on a supercritical Gaus- 
sian barrier where U = U c = 3.26m for ma = 1 (cf. 
U c = m + ^/m 2 + tt 2 /4a 2 — 2.86m for ma = 1 for a 
square barrier). Again there is a transmission resonance 
demonstrating complete penetration of the barrier. 
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FIG. 1. The zero momentum wave function for the square 
barrier V = U c (x), depicted by the heavy line. The solid 
line is the upper component and the dashed line is the lower 
component. 




FIG. 2. The zero momentum wave function for the Gaus- 
sian barrier V = U c exp(—m 2 x 2 ), depicted by the heavy line. 
The solid line is the upper component and the dashed line is 
the lower component 



While the wave functions in this case have a similar 
form to those for the square barrier, note the two turn- 
ing points which occur in the top component of the Gaus- 
sian wave function. These correspond to the points ±Xk 
where V(xk) = E + m = 2m at zero momentum. Hole 
states can propagate under the potential without expo- 
nential suppression from — xk to +xk thus demonstrat- 
ing Klein tunnelling. Note also that the condition for 
hole states to propagate under a potential of finite range 
is V > 2m which will in general not be sufficient for su- 
percriticality (we have seen for a square barrier of range 
a with ma — 1 that V c — 2.86m). So Klein tunnelling 
should exist even for subcritical potentials as was pointed 
out by Jensen et al pljj . 
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